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Abstract 

We have developed a detailed Monte Carlo simulation program for the Belle TOF 
system. Based on GEANT simulation, it takes account of all physics processes in the 
TOF scintillation counters and readout electronics. The simulation reproduces very 
well the performance of the Belle TOF system, including the dE/dx response, the 
time walk effect, the time resolution, and the hit efficiency due to beam background. 
In this report, we will describe the Belle TOF simulation program in detail. 
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1 Introduction 



A time of flight (TOF) detector[l,2] system using plastic scintillation counters 
is very powerful for particle identification (PID) in e + e~ collider detectors [2,3]. 
A TOF system with 100 ps time resolution is effective for particle momenta 
below about 1.2 GeV/c, which encompasses 90% of the particles produced in 
T(4S') decays. It can provide clean and efficient flavor tagging for CP violation 
studies in B meson decays in a B-factory experiment. 

A reliable Monte Carlo (MC) simulator is a good tool to understand the 
detector system and is important for minimizing experimental systematic un- 
certainties for physics analysis. Although one simple simulation method is a 
parametrization based on real data, it doesn't provide precise simulation for 
complex experimental situations, such as multi-track hits, hadronic interac- 
tions in the TOF counter, and shower leakage from the calorimeter outside of 
the TOF system. 

We have developed a full detector simulator for the TOF system. Based on 
GEANT3 [4] , it includes the simulation of all physics processes involving scin- 
tillation light production and its propagation in the TOF counters, the pro- 
duction of signal pulses in the photomultiplier tubes (PMT), the electronics 
process in the discriminator and the charge to time converter (QTC), and the 
effects of beam background. 

In this report, we describe the Belle TOF detector simulation. In section 2, 
we briefly introduce the Belle TOF system. In section 3 we describe the 
simulation model and the main algorithm for each process. Finally, we compare 
the simulation results with experimental data in section 4. Section 5 contains 
conclusions. 



2 TOF system in the Belle detector 

The Belle detector is a large solid- angle spectrometer based on a 1.5 T super- 
conducting solenoid magnet; a detailed description can be found elsewhere [2]. 
Figure 1 shows the Belle detector. 

The TOF system consists of 64 TOF/TSC modules located at a radius of 
1.2 m from the interaction point, which cover a polar angle of 33.7° ~ 120.8°. 
The design of one TOF module is shown in Figure 2. A module consists of 
two TOF counters of 4 cm thickness, viewed at both ends by PMTs, and one 
0.5 cm thick TSC counter, viewed only at the backward (i.e. large angle) end. 
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In order to get sufficient gain in a 1.5 T field, a fine mesh photomultiplier (FM- 
PMT) [5] is used. This tube has 24 fine mesh dynodes of 2000 mesh/inch, and 
its size is 8 cm in length and 5.08 cm in diameter. Its bialkali photocathode 
covers 50% of the end area of each TOF counter. The PMT gain is adjusted to 
be about 1 x 10 6 , which requires a high voltage of about 2100 V in the 1.5 T 
field. The FM-PMTs are attached to the TOF counter ends with an air gap of 
~ 0.1 mm that selectively passes early arrival photons and reduces the PMT 
anode current. For the TSC counters, the tubes are glued to light guides that 
are attached to the sintillator. Table 1 lists the characteristic parameters of 
the TOF scintillator and the FM-PMT used in this simulation. 

Figure 3 shows a block diagram of the TOF front end electronics (TOFFEE). 
Each PMT signal is split in two. One half is sent to a charge-to-time-converter 
(QTC) for charge measurement. The other half goes to generate two output 
signals by application of different threshold levels. The high level (HL) signal 
is used to provide a gate to the QTC, and the low level (LL) signal is fed to 
a time stretcher (TS) for time measurement. 



3 Simulation procedure 

Figure 4 shows the outline of the Belle TOF simulation program. It consists 
of the GEANT part, the main TOF simulation part including the scintillator, 
the PMT and discriminator processes, and the beam background merge. 

GEANT provides the hit information for a particle passing through the scin- 
tillator, including the four vector, hit position, time, and energy loss. 

Using the hit information from GEANT, scintillation light is produced. The 
light is then propagated in the TOF scintillator until it arrives at a PMT 
and makes a signal pulse. Finally, a pulse larger than the HL threshold makes 
a gate to measure the charge and timing of the signal. The main detector 
characteristics, such as the time resolution, are simulated in this part. 

In order to take into account the effect of beam background, real data taken 
with random triggers, which represent the beam background, are merged into 
the MC output. The inefficiency due to the dead time of the QTC is simulated 
here. 

3.1 Geant simulation up to the TOF counters 

The geometrical and material configuration of the Belle detector is imple- 
mented in the Belle GEANT full simulator [2] . Not only all sub-detector com- 
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ponents, but also materials for the support, cooling structures, and readout 
cables are included. These allow realistic simulation of a primary particle that 
goes through the inner sub-detector (and materials) up to the TOF scintillator 
and of secondary particles that are produced in the inner and outer materials. 

The material of the TOF/TSC scintillator is defined as a mixture of hydrogen 
and carbon. The atomic ratio is assumed to be 1:1. The density of this material 
is 1.032. 

The GEANT hit information for a particle traversing the scintillator volume 
is recorded at discrete steps along the particle trajectory. For a particle with 
a given energy, the step size is automatically calculated by GEANT depend- 
ing on particle species and on the characteristics of the current medium [4]. 
Figure 5 shows the time at the first step and energy loss distributions for 
2.0 GeV/c ji~ which enter perpendicular to the TOF counter. 



3.2 Scintillation process 



Using the step information from GEANT, we reconstruct a continuous path 
connecting each step in the scintillator volume. Scintillation light is uniformly 
emitted along this path. The time t tra j on this path for each light emission 
event is calculated by using the four vector of the particle. 

The number of photons is proportional to the energy loss of the particle, at a 
rate assumed to be one photon per 100 eV [6]. The angular distribution of the 
emitted light is isotropic. The time profile of scintillation emission is simulated 
by using the emission time (t emit ) probability function [7], 

E(t emit ) = — — + - e ~'— t/T3 (1) 

V ermt, 1 + R \ ^ _ ^ ^ I K J 



where T\ and r 2 are two fast decay time constants, T3 is the decay time of the 
slow component, and R is the ratio of the slow to fast components. Since the 
signal timing is determined at the leading edge in the discriminator, t\ and r 2 
are the most important parameters for TOF simulation. We use t\ = 0.9 ns 
and r 2 = 2.1 ns; these values are chosen so that this function has the same 
rise time as the BC408 scintillator light output. The values r 3 = 14.2 ns and 
R = 0.27 are given in reference [7]. Figure 6 shows the E(t emit ) distribution. 
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3.3 Light Propagation 



The scintillation light is propagated by total internal reflection and arrives at 
the PMTs located at the ends of the scintillator. We assume that light which 
has an emission angle within the trapping angle (9 trap ) can traverse the air 
gap to reach the PMT, while the remaining light escapes from the scintillator 
volume without total reflection or is lost during the endless total reflection 
between the two ends of the TOF scintillator (see Figure 7-(l)). The angle 
6 trap is given by 

trap = sin' 1 , (2) 

\Tl scint ' 



where n air and n scint are the refractive indices of the air and scintillator, re- 
spectively. 

By assuming perfect internal reflection on each of the four surfaces of the 
TOF counter, the light propagation length in the counter is calculated (see 
Figure 7- (2)) as 



Ipro cosQ i (3) 



where d is the distance from the light emission position to the end surface of 
the scintillator and 9 is the polar angle of the emitted light. This technique 
effectively reduces the CPU time, while the result is equivalent to a full sim- 
ulation of light propagation. Figure 8 shows the path length distribution of 
light which is produced at z = and reaches the forward PMT. 

The light propagation time in the scintillator is obtained by 

tpro j (4) 

V scint 

where v scint = c/n scint is the velocity of light in the scintillator. 
Attenuation in the TOF counter is simulated according to 

R(l pro ) = e~ lpro/Xpro (5) 



where \ pro is the attenuation length in the scintillator material. We use \ pro = 
400 cm, which is measured with real data. 
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The mismatch of the rectangular cross section of the TOF counter to the round 
PMT results in light loss (see Figure 7-(3)). In the simulation, the number of 
photons is reduced by the ratio of the TOF end area to the effective area of 
the PMT photocathode. 



3.4 PMT response 

We consider the angular response (R an gie) of PMT sensitivity, which is a de- 
pendence on the angle of incident light on the photocathode [8] , 



where 6i is the light incident angle at the PMT front glass. We find 0.66 for 
the efficiency averaged over all arrival directions at the PMT. The number of 
photons is reduced by this factor. 

At the PMT photocathode, photoelectrons are produced with a probability 
which is called the quantum efficiency The average quantum efficiency of 
the photocathode is calculated to be 0.21 for the spectrum from the BC408 
scintillator. Also, the number of photoelectrons is reduced by a factor of 0.6, 
called the collection factor, due to the fine mesh dynode structure. We reduce 
the number of generated photoelectrons by these factors. 

The transit time of a single photoelectron signal in a FM-PMT, t TT , is mod- 
eled by a transit time(TT) of 5.5 ns and a r.m.s. transit time spread (TTS) of 
0.31 ns, which are measured for the R6680 FM-PMT. By adding the smeared 
transit time to the photon timing information described in the previous sec- 
tion, the pivot time of single photoelectron signal at the anode is obtained, 



Figure 9-(l) shows the time distribution of the photoelectrons at the anode 
for a typical event, a 2.0 GeV/c fi~ incident on a TOF counter at z — 0. 

To simulate the time walk effect due to the leading edge discrimination scheme 
used for the Belle TOF, a time response function Vi(t) for a single photoelec- 
tron is introduced. Using this response function, the PMT output signal is 
obtained by summing all the single photoelectron signals in the PMT for an 

event, 



Rangie COsOi 



(6) 




(7) 



VpMT{t) = £>(<), 



(8) 



i=l 
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where 



Vm = GC ° Jt*e-Wdf 



(9) 



Here G is the gain of PMT, C e is the conversion factor for charge to voltage, 
and r is the time constant. We use r = 6 ns so that the rise time of this 
function is the same as that of an R6680 PMT, and 25 ps time bin size which 
is the same as the least count of the Belle TOF readout system [1]. 

Figure 9- (2) shows a typical PMT output signal for the same event as shown 
in Figure 9-(l). The simulated signal is consistent with actual signals observed 
in the Belle TOF counters in both pulse height and rise time. 

3. 5 Simulation of readout electronics 

Each PMT signal is examined by readout electronics using double thresholds. 
A signal larger than the high level (HL) threshold provides a gate to measure 
the signal timing (TDC) which is given at the moment when the PMT pulse 
crosses the low level (LL) threshold in the discriminator simulation. We use 
the same LL threshold as the real experiment (50 mV). A HL value 50% 
higher than that in the real experiment (150 mV) is used for fine tuning of the 
hit inefficiency, which is caused by PMT-to-PMT gain variation in the real 
experiment. 

Finally, we introduce a timing uncertainty, parametrized by a Gaussian with 
a width of 45 ps, which is due to the uncertainty of the beam crossing time 
(t = 0) determination, the time jitter caused by electronic noise in the discrim- 
inators, and the time jitter in the TS readout. Figure 10-(1) shows the TDC 
distribution for perpendicular incidence on a TOF counter, and Figure 10- 
(2) shows TDC values as a function of particle hit position (z-hit) along the 
counter. The z-hit value is obtained by extrapolating the particle's trajectory 
using tracking information from the central drift chamber (CDC) [2]. The 
filled circles in Figure 10 are data and the open circles are MO There is good 
agreement. 

The charge output (ADC) is given by the time integral of the Vpmt(^) divided 
by the impedance of the PMT circuit. The width of the simulated charge distri- 
bution is narrower than that of the real data. Possible sources of the difference 
include statistical gain variations of the PMTs [9] (the simulation assumes a 
constant gain for a single photoelectron), gain variation among individual 
PMTs, and the imperfection of the simplified model of light propagation. To 
better represent the data, we adjust the charge distribution by introducing a 
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scaled Poisson distribution, 

ADC = Q orig ( Sw i N M-n pe ) + \ ? (lQ) 

where ADC is the corrected charge, S w is the width scale parameter, N{n pe ) is 
a random number given by a Poisson distribution generator which has mean 
n pe and Qorig is the original charge from the VpMrit) integration. In a com- 
parison between we found that a value S w = 2 gives the 
best agreement. Figure ll-(l) shows the resulting charge distribution com- 
pared with the real data. Figure 1 1-(2) shows the mean ADC value as a 
function of the z-h.it position. All points are in good agreement. By fitting 
\og(ADC f / ADC b) as a function of z-h.it (see Figure ll-(3)) to a linear func- 
tion, we obtain good agreement for the effective light attenuation length (A) 
along the z direction [10]. 

One of the main parameters in the TOF calibration and reconstruction is a 
time walk correction coefficient that is used to remove a pulse height depen- 
dence [1] of the digitized time. It is the coefficient of a 1/V A DC term, and is 
sensitive to the rise time of the PMT signal. Figure 12 shows the (uncorrected) 
TDC versus 1/V ADC for data and MC. They are in good agreement. 

3. 6 Beam Background Simulation 

The effect of the beam background is simulated by merging beam background 
data into the MC signal output. The background data are taken by random 
triggers during the real experiment. The typical hit rate of a single TOF PMT 
is 15 ~ 25 kHz. Because of its greater scintillator width, the TSC counter 
hit rate is more than 2 times higher than the TOF rate. The hit rates are 
expected to increase in the higher luminosity environment in the near future. 

Figure 13 presents a flow chart of the beam background merging. For each 
PMT, MC and background are merged on an event-by-event basis. First, we 
consider the overlapping of two PMT signals. According to the algorithm in 
the readout electronics, if two signals are closer than the ADC gate width 
(100 ns), we take only the earlier signal TDC and the charge of the later 
signal is added to the earlier one. Otherwise, both signals are recorded. Then 
the QTC dead time and the number of edges are taken into account. 

The major effect of the beam background is to increase the inefficiency caused 
by the QTC dead time, due to signal conversion and recovery times. It is 
simulated by discarding the later TDC signal in the case that the MC and 
beam background times are closer than the QTC dead time, which is found 
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to be 0.55 ns per ADC count (Each ADC count corresponds to a charge of 
0.25 pC.) in real data. 

Another source of inefficiency is the limit on the number of records in the 
TDC1877 system [2]. The limit is 8 edge records for each PMT in an event. 
Because the TS system produces 4 edges for one PMT signal [1], a maximum 
of 2 PMT signals can be recorded for an event. In our simulation, we count the 
number of edges for the individual TDC channels, and remove the signal(s) 
after the 8th edge in the merged data. Since we don't use the TS for the TSC 
counters, this inefficiency is not present for them. 

Figure 14 shows the TOF and TSC hit inefficiencies in /z-pair events (e + e~ — > 
as a function of the TOF trigger rate, which is defined as the rate for at 
least one TOF-TSC coincidence among the 64 TSC counters. The TOF PMT 
singles rate is about 25% of this TOF trigger rate, and the TSC singles rate 
is about two times the TOF singles rate. Here the hit inefficiency is defined 
as the fraction of muons predicted to hit a TOF counter for which there is no 
associated counter signal. In the figure, the solid circles are for data and the 
open circles are for MC with beam background. The hatched part indicates 
the inefficiency given by the original MC without beam background. This 
intrinsic inefficiency is due to the dead space between counters. The higher 
TSC inefficiency is due to the higher background rate. The hit inefficiencies 
for MC and data agree to within 0.3%. 



4 Result 

Figure 15 shows the dt = T^ rect — T pred distribution for all available TOF hits 
in a /i-pair event sample. Here T pre d is the time of flight predicted using the 
track path length calculated from the fit to hits in the CDC, and T^"' ect is the 
weighted average of the observed times in the forward and backward PMTs 
after correction for the z-hit position and time walk [1]. The distribution is 
fitted to a symmetric Gaussian to determine the time resolution. The overall 
time resolution (100 ps) is well represented by the MC. Figure 16 shows the 
time resolutions for forward and backward PMTs and for the weighted average 
time as a function of z-hit . 

The PID from the Belle TOF system is based on a ratio of likelihoods for two 
particle hypotheses. For Kj n separation, a particle is identified as a kaon using 
the likelihood ratio PID(K) = P K j{P^ + P K ), where the P K>W are the PID 
likelihoods for the particle hypotheses in the subscripts. The PID performance 
is evaluated using the charm decay D* + — > D tt + , followed by D° — > K~tt + . 
Figure 17 shows the measured K ID efficiency and 7r fake rate as a function 
of momentum, where PID(K) > 0.6 is required to identify a kaon. The MC 
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reproduces the K ID efficiency and fake rate as functions of momentum with 
maximum deviations of 0.7% and 1.8% , respectively. 



5 Conclusion 

In this report we describe the full detector simulation for the Belle TOF sys- 
tem. We developed a GEANT-based simulator for the processes in the scintil- 
lator and the PMT, the readout electronics of the discriminator and the QTC, 
and the effect of beam background. 

Scintillation light is produced using the hit information from GEANT. Using 
simple models for light propagation in the TOF scintillator and for the PMT 
response, we accurately reproduce the real PMT signals. The simulation repre- 
sents the detector characteristics of the TOF system such as light attenuation, 
the ADC distribution, the effective light velocity, the TDC distributions and 
the time resolution. Through a beam background merging process, the QTC 
dead time and the limit on the number of records in the TDC1877 system 
is simulated, so that the TOF hit inefficiency dependence on the background 
rate is accurately modeled. Since the dependence on pulse height and on z-hit 
is well reproduced, the simulation is a good tool for the investigation and de- 
velopment of calibration and reconstruction methods. The K ID efficiency and 
its momentum dependence are in good agreement with the data. The TOF 
simulation provides an important component of the Belle MC used for physics 
analysis. 
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Table 1 

Characteristic parameters of the TOF scintillator and the FM-PMT. 



TOF scintillator (BC408) 



Base 


Polyvinyltoluene 


Density 


1.032 


Refractive index 


1.58 


Rise Time 


0.9 ns 


Decay Time 


2.1 ns 


Pulse Width 


~ 2.5 ns 


Light Attenuation Length 


~ 300 cm 


Wavelength of Maximum Emission 


425 nm 


FM PMT (R6680) 


Effective photocathode diameter 


39 mm 


Transit Time Spread 


320 ps (r.m.s) 


Quantum Efficiency 


~ 0.21 


Electron Collection Factor 


0.6 


Rise Time 


3.5 ns 


Fall Time 


4.5 ns 


Pulse width 


6 ns (FWHM) 
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Fig. 1. A cut-away view of the Belle detector. 
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Fig. 2. A TOF/TSC module. 
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Fig. 3. A block diagram of the TOF front end electronics (TOFFEE). 
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Fig. 4. A block diagram of the simulation program. 
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Fig. 5. The first-step time (left panel) and energy loss (right panel) distributions. 




Fig. 6. The emission time probability for scintillation light. The line is the function, 
the histogram is the simulated result. 
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7. The light propagation in a scintillator counter. 
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Fig. 8. The light propagation length distribution. 
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Fig. 9. (1) The time distribution of the photoelectrons, (2) the PMT pulse output 
(with negative sign), and (3) the PMT response for a single photoelectron. 
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(2) TDC vs. z-hit 
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Fig. 10. (1) The TDC distribution for -5< z-hit<5 cm, (2) the TDC distribution as 
a function of z-hit. The solid circles are data; the histogram and open circles are 
MC. 
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Fig. 11. (1) The ADC output distribution. (2) The mean ADC value as a function 
of z-hit. (3) The dependence of log(ADCp / ADCb) on z-hit. The solid circles are 
data; the histogram and open circles are MC. 
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Fig. 12. The variation of the TDC output vs 1/VADC for -5< z-hit<5 cm. The 
solid circles are data and the open circles are MC. The TDC and ADC values are 
for the forward PMTs. 
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13. A flow chart for the merger of beam background with simulated events. 
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Fig. 14. The hit inefficiency as a function of the TOF rate. The solid circles are data 
and the open circles are MC. 
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Fig. 15. The dt distribution. The circles are data; the histogram is MC. 
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Fig. 16. The TOF resolution as a function of z-hit. The top panel is MC; the bottom 
is data. 
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Fig. 17. The K ID efficiency and ir fake rate as a function of momentum. The circles 
are data; the histogram is MC. 
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